Griineisen ratio at the Kondo breakdown quantum critical point 
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We show that the scenario of multi-scale Kondo breakdown quantum critical point (QCP) gives 
rise to a divergent Griineisen ratio with an anomalous exponent 0.7. In particular, we fit the 
experimental data of YbRh2{Sio.95Geo.o5)2 for specific heat, thermal expansion, and Griineisen 
ratio based on our simple analytic expressions. A reasonable agreement between the experiment 
and theory is found for the temperature range between 0.4 K and 10 K. We discuss how the Griineisen 
ratio is a key measurement to discriminate between the Kondo breakdown and spin-density wave 
theories. 
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Heavy-fermion quantum criticality is a typical example 
of a quantum system where both strongcorrelations and 
Fermi surface effects play a major role [l| . The standard 
model of quantum criticality in a metallic system is a z = 
2 critical theory, often referred as Hertz-Moriya-Millis 
theory Q], where z is the dynamical exponent relating 
the variation of the energy with the momentum, lo ~ 

. Unfortunately, many heavy fermion compounds have 
been shown not to follow the spin-density-wave (SDW) 
theoretical framework 0, |3) IE H ■ 

An interesting suggestion is that the heavy-fermion 
quantum transition is analogous to the Mott transition 
0, H, • The arguments in support of this view are 

the divergence of the effective mass near the QCP [sj and 
the presence of localized magnetic moments at the tran- 
sition towards magnetism [J]. Combined with the Fermi 
surface reconstruction at the QCP ^6'], this quantum 
transition is assumed to show a breakdown of the Kondo 
effect as an orbital selective Mott transition, where only 
the f-electrons experience the metal-insulator transition. 

Recently, thisproblem has been re- visited in the slave- 
boson context [3, . The main idea is that the Kondo 
breakdown QCP is multi-scale. Dynamics of hybridiza- 
tion (holon) fluctuations is governed by spinon-electron 
polarization. An important observation is that there 
should exist a Fermi surface mismatch q* = \kp — kp\ 
between Fermi momentum kp for spinous and fc|, for 
conduction electrons since fillings of spinous and elec- 
trons differ from each other. Fermi surface mismatch 
gives rise to an energy gap E* for such spinon-electron 
fluctuations. Although it depends on the value of g*, this 
energy scale is shown to vary from ©(lO*^) mK to O(IO^) 
mK. When E < E*, holon fluctuations are undamped 
without considering gauge fluctuations, thus described by 
z = 2 dilute Bose gas model. On the other hand, when 
E > E* , holon fluctuations are dissipative since spinon- 
electron excitations are Landau damped, thus described 
by z = 3 critical theory. Based on the z — S quantum 
criticality, recent studies [§, [l3| have found quasi-linear 
electrical transport and logarithmically divergent specific 
heat coefficient in d = 3, consistent with an experiment 

i- 

In this paper we study the Griineisen ratio (GR) 



rs(r, T) = as{r,T)/cs{r,T) based on the multi-scale 
Kondo breakdown QCP scenario where as{r, T) — 

^ drdT'^ and Cs(r, T) = -T^-^^^ are thermal ex- 
pansion and molar specific heat with molar volume v, 
respectively. Recently, GR has been proposed as one pos- 
sible measure characterizing the nature of a QCP [ll| . A 
remarkable feature is that the GR diverges at any QCP 
with an anomalous exponent depending on the nature of 
the quantum transition. Consider the scaling expression 
fs{r,T) = b-'^'^+^^frirb^/" ,W) for the free energy near 
a QCP at spacial dimension d, where r k. ^p^" is a dis- 
tance from the QCP (Pc) with a pressure-unit constant 
Pq and 6 is a scaling parameter with a correlation-length 
exponent ly. Evaluating the thermal expansion and mo- 
lar specific heat, we find the following scaling expressions 



with Gr — i^{d — uqz) /uq where i/q is an exponent associ- 
ated with the third law of thermodynamics , and 



with Gt 



The GR at the 



Ts{T;r = 0) = -GT T"^ 

1 [dfAt,l)/dt]t=a 
Pqv vd{d+z) /r(0,l) 

QCP, Ts{T]r — 0) exhibits the scaling exponent x — ^ 
in any dimension. 

In a recent measurement 12] it was reported that 
for 1^6i?ft.2(<S'jo. 95^60.05)2, the specific heat coefficient 
can be well fitted by 7(r) = C{T)/T cx In(^) for 
- Q.'iK < T < IQK, with an energy scale w "iOK 
identified with its Kondo temperature, and 7(T) oc T^'^/^ 
for T <~ 0.3i^. The thermal expansion coefficient was 

fitted as a{T)/T oc - In(^) for IK < T < IQK with a 

temperature scale ~ and a(T)/T cx ao + ai/T 

for ^ OAK < T < IK with ao « 3.4 x IQ-^K-'^ and 
ai ~ 1.34 X lO^^A'^^. Finally, the experiment shows 
that the GR diverges with an exponent x ^ 0.7 ± 0.1. 
This invalidates the SDW scenario, since we have x = 1 
owing to z = 2 and v — 1/2, where this critical theory is 
beyond its upper critical dimension in d > 2. 

In this study we show that the scenario of multi-scale 
Kondo breakdown QCP gives rise to a divergent GR with 
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SDW Kondo breakdown 



FIG. 1: (Color online) Schematic diagram of Griineisen ratio 
Ts{r; T ^ 0) for both the SDW and Kondo breakdown scenar- 
ios, where AF and PM represent antiferromagnetic and para- 
magnetic phases in the SDW context, and SL and HF denote 
spin liquid and heavy fermion states in the Kondo breakdown 
scenario. Two kinds of boson excitations, hybridization (Fg) 
and gauge (Ff ) fluctuations, contribute to Griineisen ratio in 
the Kondo breakdown scenario. Gauge fluctuations play an 
important role for thermal expansion in the HF phase more 
than in the SL phase. As a result, asymmetry is expected 
to appear for Griineisen ratio around the Kondo breakdown 
QCP. 



the exponent 0.7. In particular, we fit the experimental 
data [l^ of F6i?/i2(5'io. 95660.05)2 for specific heat, ther- 
mal expansion, and GR with simple analytic expressions 
[Eqs. (3) and (4)], for which the asymptotic behavior is 
summarized in Table I. The z = 3 quantum criticality in 
d = 3 turns out to play an essential role for thermody- 
namics near the QCP of F5i?ft,2('S'io. 95660.05)2- 

a,{T) csjT) r,(r) 

TABLE I: Thermodynamics in the z — 3 regime {d = 3) 

Interestingly, varying an external parameter gives an 
opportunity to distinguish the Kondo breakdown from 
the SDW scenario (see Fig. 1). It has been shown 
that the thermal expansion coefficient should change sign 
across the SDW QCP in the zero temperature limit, and 
GR also does accordingly In the Kondo breakdown 
scenario two kinds of collective excitations, hybridiza- 
tion and gauge fluctuations, contribute to thermal ex- 
pansion. Hybridization fluctuations give rise to the same 
sign change as the SDW fluctuations while gauge fluctua- 
tions do not. Considering that gauge fluctuations should 
remain gapless in the spin liquid phase due to gauge in- 
variance, their contribution for thermal expansion is van- 
ishingly small in the spin liquid phase. On the other 
hand, they contribute to thermal expansion heavily, ap- 
proaching the QCP in the heavy-fermion phase owing 
to the Anderson-Higgs mechanism. Taking into account 
both hybridization and gauge fluctuations, an asymmet- 
ric feature of GR is expected to appear around the Kondo 



breakdown QCP. 

We start from the U(l) slave-boson representation of 
the Anderson lattice model in the large-t/ limit 

^ = Yj ^^'^(^r - Ai)ei<x -tJ2i4a('3<y + H.C.) 

+Vj2ibjlc^a + H.C.) + h\drh 

i i 

+ E fl^^r + 6/)/.. + ^ Y.^flXvf,. + H.C.) 

i {ij) 

+1 J2 + flu - 1) + ^ E i^^^i'- (1) 

« (ij) 

Here, Ci^ and d^o- = ^j/io- are conduction electron with 
a chemical potential n and localized electron with an 
energy level e/ respectively, where bi and are holon 
and spinon, associated with hybridization and spin fluc- 
tuations. The spin-exchange term for the localized or- 
bital is introduced for competition with the hybridization 
term, and decomposed via exchange hopping processes of 
spinons, where Xij is a hopping parameter for the decom- 
position. Xi is a Lagrange multiplier field to impose the 
single occupancy constraint 6-6,i -I- fj^fia = ^/2, where 
N is the number of fermion flavors with a = \, N . 

The slave-boson mean-field analysis has shown an 
orbital selective Mott transition [lj| as breakdown of 

Kondo effect at J « Tr-, where Tk = D exp ^ j^p^y2 ^ 

is the Kondo temperature with the density of states Pc 
for conduction electrons [9]. If we try to understand the 
GR in this level of approximation, we find x = 1 for the 
GR exponent. Actually, one can check that the mean- 
field free energy satisfies the following scaling behavior 
fMF{B, T) = T^J'iBT-^), where T{x) is an analytic 
function and B — V{b) is the effective hybridization. 

Fluctuation-corrections are important at the Kondo 
breakdown QCP, where both hybridization and gauge 
fiuctuations should be taken into account carefully. Such 
fiuctuations are treated on an equal footing in the Eliash- 
berg framework, where momentum dependence in self- 
energies and vertex corrections are neglected, justified by 
the Migdal theorem and large N approximation [l^ . 

For a systematic study of thermodynamics, we con- 
struct a Luttinger-Ward (LW) functional in the Eliash- 
berg framework, composed of contributions from conduc- 
tion electrons, spinons, holons, gauge fluctuations and 
their self-energy parts. One can derive self-consistent 
Eliashberg equations for the self-energies from variation 
of the LW functional with respect to each self-energy. 
Using these equations, one is allowed to simplify the LW 
functional as Flw = FpL + ^fl + Pb + Fa, where the 
first two parts represent Fermi liquid contributions for 
conduction electrons and spinons while the latter two 
parts express hybridization and gauge contributions, re- 
spectively. Such fermion contributions are sub-dominant 
compared with boson contributions, and they can be ig- 
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nored in the low energy limit. Accordingly, thermal ex- 
pansion and specific heat can be approximated as follows 
near the Kondo breakdown QCP, a,(T) « ab(T) + aa(T) 
and Cs(T) « Cb{T) + Ca{T), respectively. As a result, the 
GR is found to be r,(r) « £^g±^^. 

The bosonic part of the free energy is given by 



Here, 7^ = ^ and 7^ = ^ + ^""^ ''"^'^ are Landau 
damping coefficients for holon and gauge fluctuations, re- 
spectively, where Vp is the Fermi velocity for conduction 
electrons, a = is an effective ratio between band- 
width of each fermion sector, and fd is associated with an 
ultra-violet cutoff for gauge fluctuations. A;, is the mass 
for the hybridization fluctuations, identifying the Kondo 
breakdown QCP with Af, = 0. Aa is the mass for gauge 
fluctuations, resulting from Anderson-Higgs mechanism, 

thus related with the mass of holon as A,, = ^^^^^ At, 
where W{, is the strength of local interactions for holons, 
phenomenologically introduced. Fc is the condensation 
part. 

Performing the frequency summation and momentum 
integral, we find the specific heat and thermal expansion 
coefficients at the QCP, 



c,iT>E*) 
T 

Cs{T < E*) 
T 



and 



T 

asiT < E*) 



— c 



dP 
dAb 
dP 



2^1 E*-^ + 



(4) 



where Cr 



ycothy 



and Ca — 

'^y^i^) are positive numerical 
constants. The condensation part is assumed to be al- 
most constant for temperature dependence, thus can be 
ignored for thermal expansion. Note that there is an 
with pressure P in the thermal 



unknown constant 



OA, 
dP 



expansion coefficient, determining its overall sign. Re- 
calling that it is negative for F6i?/i2(5'jo. 95^60.05)2 [l^ - 
we see < 0. This implies that pressure puts the 
QCP of y5i?/i2(S'«o. 95^60. 05)2 toward the heavy-fermion 
side if it is identified with the Kondo breakdown QCP. 
One can check that holon thermodynamics is consistent 
with z = 3 scaling for T > E* while gauge thermody- 
namics is for all temperatures. In addition, both specific 
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FIG. 2: (Color online) Specific heat coefficient, where the 
blue dotted line represents an experimental data and the red 
thick line does our theory. Inset: Specific heat contributions 
from hybridization fluctuations (red) and gauge fluctuations 
(green). Note that the speciflc heat coeflicient from hybridiza- 
tion fluctuations is twice bigger than that from gauge fluctu- 
ations. 




FIG. 3: (Color online) Thermal expansion coefficient, where 
the blue dotted line represents an experimental data and the 
green thick line does our theory. Inset: Thermal expan- 
sion contributions from hybridization fluctuations (black) and 
gauge fluctuations (green). Note that the thermal expansion 
coefficient from gauge fluctuations is much larger than that 
from hybridization fluctuations although this physics depends 
on the local-interaction strength ut for hybridization fluctua- 
tions. 



heat and thermal expansion coefhcients are constant for 
holon fluctuations at T < E*, consistent with Fermi liq- 
uid physics. 

Using Eqs. (3) and (4), we try to fit the experimen- 
tal data of Ref . ■ The density of states Pc — jjy for 
conduction electron with the bandwidth D = 6t ^ W^K 
and the ratio between bandwidths a — 10^^ are fixed. 
Considering Vp = ^ « 2t for the sphere Fermi surface, 
7h is determined as 7b = 2.693. V can be deduced from 
the Kondo temperature Tk ~ •iOK with ej — — -j; thus 
V = 0.293D. fd is used as a fitting parameter, determin- 
ing "fa — 1.367. The cutoff A in specific heat is set to be 
an effective bandwidth for localized spins, i.e., A « aD. 
E* is approximately given by the upturn temperature for 
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FIG. 4; (Color online) Griineisen ratio, where the blue dot- 
ted line represents an experimental data and the red thick line 
does our theory. Inset: Comparison between the d = 3 Kondo 
breakdown (red) and d = 2 SDW (black) theories, where the 
SDW theory exhibits more singular behavior at low temper- 
atures. 

specific heat, here E* « O.SfT. For thermal expansion, 
we have two free parameters, Uh and ^f- 

Fig. 2 shows the fitting for the specific heat coefficient 
Cs{T)/T. For T > E* , we have a very good matching 
unlike for T < E* . Although the origin of this upturn 
behavior is not explained, yet clearly, two dimensional 
ferromagnetic fluctuations [14| may give one possible ex- 
planation, resulting in cjrM{T)/T oc T~^^^. From the in- 
set figure, we can conclude that both hybridization and 
gauge fluctuations are important for specific heat near 
the Kondo breakdown QCP. Considering that the hy- 
bridization fluctuations arise from collective excitations 
of conduction electrons and spinous, and gauge fluctua- 
tions result from spinon current-current correlations, one 
can expect that both fluctuations will contribute in a 
similar fashion. 

Fig. 3 shows the fitting for the thermal expansion co- 
efficient as{T)/T, where we have a rather good agree- 
ment between experiment and theory. Although we 
have used two free parameters Ub and ^r, such pa- 



rameters can change only the overall scale, thus one 
may regard that only one parameter is used. The in- 
set figure exhibits that contributions from gauge fluctua- 
tions are much larger than those from hybridization ones. 
Although this physics depends on the local-interaction 
strength Uh for hybridization fluctuations, it is valid as 
far as Ub T^'^^' satisfled in Eq. (4), i.e., in the 
weak coupling limit preserving the present picture of the 
Kondo breakdown QCP. Our fitting for the thermal ex- 
pansion coefficient may be the first explicit demonstra- 
tion, supporting importance of gauge fluctuations. 

Fig. 4 shows a reasonable match between experiment 
and theory for the GR above the upturn temperature. As 
shown in the inset figure, the d = 2 SDW theory shows 
more singular behavior at low temperatures, deviating 
from the experiment more seriously. 

In conclusion, we have fitted the experimental data of 
y6i?/i2 (5*10.95060.05)2 for specific heat, thermal expan- 
sion and Griineisen ratio based on simple analytic formu- 
lae in the multi-scale Kondo breakdown scenario. Both 
hybridization and gauge fluctuations contribute to spe- 
cific heat in a similar fashion around the QCP. Gauge 
fluctuations are more important in the heavy-fermion 
phase than in the spin liquid phase for thermal expan- 
sion, causing an asymmetry for Griineisen ratio around 
the QCP. This feature can be used to discriminating the 
Kondo breakdown scenario from the SDW framework. 
These z — 3 critical fluctuations explain the divergent 
Griineisen ratio with the anomalous exponent 0.7 beyond 
the SDW theory. We suggest that two dimensional ferro- 
magnetic fluctuations may give one possible explanation 
for thermodynamics in the low temperature region below 
E* , not captured in the present framework. 
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